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^ Abstract 

^ We study the decay B -^ K*£^£~ in the QCD factorization approach and propose 

^ a new integrated observable whose dependence on the form factors is almost negligible, 

K> consequently the non-perturbative error is significantly reduced and indeed its overall 

H theoretical error is dominated by perturbative scale uncertainties. The new observable we 

propose is the ratio between the integrated forward-backward asymmetry (FBA) in the 
[4,6] GeV^ and [1,4] GeV^ dilepton invariant mass bins. This new observable is particu- 
larly interesting because, when compared to the location of the zero of the FBA spectrum, 
it is experimentally easier to measure and its theoretical uncertainties are almost as small; 
moreover it displays a very strong dependence on the phase of the Wilson coefficient Cio 
that is otherwise only accessible through complicated CP violating asymmetries. We il- 
lustrate the new physics sensitivity of this observable within the context of few extensions 
of the Standard Model (SM), namely the SM with four generations (SM4), an MSSM with 
non-vanishing source of flavor changing neutral currents in the down squark sector and a 
Z^ model with tree level flavor changing couplings. 



1 Introduction 

The importance of the flavor-changing semileptonic decays, b -^ s{d)l^l~ for searching effects of 
new physics has been recognized for a very long time [ll-lil.The process is especiaUy important 
because its amphtude is very sensitive to the presence of heavy virtual quarks, scaling as 
(massquark)^. The inclusive BR{B -^ Xsi^i') = (1.60 ± 0.51) x 10"^ is measured (5}j7| and is 
consistent with the NNLO prediction of the Standard Model (SM), (1.62±0.14) x 10"Kvith an 
accuracy of around 9% f8j[9] ; given the large experimental uncertainty, this measurement allows 
for new physics contributions at the ^ 30% level. The reaction is also of great interest as it 
oflFers numerous related observables through which to test the SM precisely and to discover new 
phenomena. The dilepton invariant mass (q^) spectrum, forward-backward asymmetry (FBA) 
as function of q'^^ the location of the zero in the FBA and the possibility of searching for CP 
violation via, for example, partial rate asymmetry (PRA) are some of the attractive features of 
this reaction. As usual while theory predictions tend to be more reliable for inclusive modes, 
experimentally the related exclusive modes such as 5 ^ K{K*)i^i~ ^ Bs -^ (l)i^i~ are more 
readily accessible. A particular challenge for theory regarding these exclusive modes flO] is their 
dependence on form factors, which are manifestly properties of bound states and therefore of 
non-pert urbative character. B factories of course have been able to study both inclusive and 
exclusive modes but at the LHCb inclusive studies seem rather diflicult. On the other hand 
many of the exclusive modes are quite distinctive and in particular B -^ i^*^+^~ , Bg -^ 0^+^" 
should be a focus of intense study at the LHCb and at the super B factories. 

In this paper, we will focus on the forward backward asymmetry in the q'^ distribution of 
B -^ K''£^£~ ^ Bs -^ (j)£^£~ . The topic has gained renewed interest recently as measurements 
at BaBar fTI|, Belle fT2| and CDF |13| all seem to show weak indications of deviations from 
expectations of the SM. Lack of statistics in these flrst observations may be the underlying 
reason but the prediction of the SM also has considerable uncertainties and that is a cause 
of some concern to which we will try to address in this work. As already mentioned the 
dominant source of the uncertainties is the form-factor dependence; these are not amenable to 
perturbation theory and improvement using non-perturbative methods tends to be rather slow. 

Here we propose the ratio of the integrated forward-backward asymmetry in the [4, 6] GeV^ 
and [1,4] GeV^ q'^ regions (just above and below the location of the zero of the spectrum). In 
this ratio, the form factors uncertainties cancel to a great extent and the total error is controlled 
by scale uncertainties whose perturbative origin allows for future improvements. 

2 B ^ K*ii in QCD factorization 

The effective Hamiltonian responsible for the short distance b -^ s£^£~ transition is 
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where the definition of the operators and the corresponding SM matching conditions can be 
found in Ref. 8p4|*. Primed operators are obtained from the unprimed ones by the substitution 



L ^ R. The electroweak penguin and O5 operators contributions to B ^ K^'M can be neglected 
because their impact is much smaUer than the overaU theory uncertainty. 

It is possible to show that, for small values of the dilepton invariant mass (g^<6 GeV^), 
the B -^ K'^iM matrix element of all the operators appearing in Eq. dl]) obey the following 



factorization formula 15 : 



{K:£i\0,\B) = CI Uq') + E /d^ d" T^a,±iQ^u,^) ^B,±i^) ^K*{u) + O (A^coM) (2) 

where a =_L, || refers to the transverse (longitudinal) polarization of the X*, ^a are soft form 
factors and $^ ^ and $^* are the light cone wave functions of the B and K* mesons. The 
quantities C^ and T^ ^ can be calculated in perturbation theory. 

B -^ K^'M matrix elements can be decomposed using a basis of eight transversity amplitudes 
(see, for instance, the discussion in Sec. 3.2 of Ref. |14|): A±l^ji^ ^\\l,Ri ^0L,i?, A and As- In 
this notation, Aj_ and A\\ correspond to the two possible transverse polarization states of the 
K* and are described by the form factor ^j_] Aq, At and As involve a longitudinally polarized 
X* and are described in terms of ^\\. 

Up to electromagnetic corrections, the various B -^ K^'M matrix elements of the semileptonic 
and magnetic moment operators O7 g ^^ql ^an be exactly expressed in terms of the seven form 
factors V{q^)^ ^1,2,3(9^) and Ti^2,3{q'^)- Up to power corrections, these form factors can be 
written as the sum of soft and hard factorizable contributions using Eq. ^. It is important to 
realize that the separation between soft and hard contributions is subject to a certain degree 
of arbitrariness; therefore it is necessary to choose a factorization scheme and adopt precise 
definitions of the soft form factors in terms of the full QCD ones. Following the analysis of 
Ref. n!6] we adopt the following definitions: 



The form factors V^ A\ and A2 are taken from the Light-Cone Sum Rule analysis of Refs. [Tip!?] . 
In the numerical analysis we vary independently the g^ = values of the two soft form factors 
^^ and (^11 over the one-sigma errors given in Table fl} The matrix elements of all the other 
operators do not admit a simple interpretation in terms of form factors but can be expressed 
using Eq. ([2]) in which ^^ and ^\\ are the same soft form factors introduced in Eq. ([3]). 



*Note that here we adopt a more conventional definition of the semileptonic operators Og^io and, unlike 
Ref. [81, we factor out Qfe/47r in their definition so that their Wilson coefficients start at order 0(a^,ag). 

■^For definiteness, we will continue to explicitly mention B -^ K^'t^t' only, however, it should be obvious 
that our equations and the ensuing discussions are applicable equally well to Bs -^ (j)£~^£~ . 

''"This applies also to the matrix elements of the scalar and pseudoscalar operators Og p. 



Let us now sketch how aU the different parts of this calculation come together. Up to 0{as) 
corrections, the three most relevant transversity amplitudes are schematically given by 

A^L,R cc [{Cf + a)T{C,, + C[,)]V{q') + ^^^{C, + C',)T,{q') (4) 

A\\L,R oc [{Cf-C',)T{C,,-C[,)]A,{q') + '^^^{C,-C',)T,{q') (5) 

AoL,R OC [iCf-C',)TiC\,-C[,)]{A,iq'),A2iq')} + iC,-a,){T,iq%T,iq')} (6) 

where {•, •} indicates a linear combination, Cf^ = Cq + Y{q'^) and Y{q'^) is the sum of the 
leading order matrix elements of the current-current and QCD penguin operators Oi-q and 
can be found for instance in Ref. |15j. All Wilson coefficients are evaluated at a scale /i^ ^ m^. 



Complete expressions for the eight transversity amplitudes can be found in Refs. [14, 18, 19 

Part of the 0{as) corrections are buried inside the form factors; therefore the first step 
consists in replacing each form factor with the corresponding QCD factorization expansion 
(e.g. V = ^± {rriB + mK*)/mB)- The complete expansion of the form factors in terms of ^j_^\\ 
can be found in Refs. fl5l|T6l[20] . In a second step all the remaining non-factorizable corrections 
originating from the remaining operators (other than Or^g^io) have to be included. Since the 
latter can always be described in terms of B ^ i^*7* matrix elements it is customary to lump 
them together in some effective photonic form factors: 

{C,-C',)Uq') ^ rr{q') (7) 

(C, + C',) T^iq') ^ r+(cz2), (8) 

where the index ± refers to the combinations C^±C^^. The quantities r- are obtained from the 
Ti introduced in Ref. fTs] with the appropriate C7 ^ C^± Cij substitution. Following Ref. [Ts] , 
the contribution from the Y[q^) term in C^^ is also included in the r^; therefore we also have 
to replace Cf^ -^ Cc 
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Note that the factorization scheme that we adopt implies the absence of 0{as) corrections 
to the form factors V{q'^)^ ^i(?^) and ^2(9^) (that are therefore simply expressed in terms of 
the soft form factors ^±^\\ without any 0{as) contribution). The form factor A^^q'^)^ on the 
other hand, receives non-trivial 0{as) corrections. Effects on the photonic form factors Ti{q^) 
are included in the definition of the r^. 

Once the transversity amplitudes have been calculated at order a^, it is trivial to express 
the fully differential B -^ i^*(^ K7t)££ decay width. In the limit m^ -^ 0, one finds: 



where 



dr 


= ^ (11 + 2/| + (2/| - ID cos^ e, + 


dg2 d cos 6i 


I2 -- 


= \(^Ai^+A\^+Al^+Af^) 


It - 


- A"P + K' 



(9) 



(10) 

(11) 
Fq = 2Re{AlA^* -AfAf*) (12) 
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mf = (4.68 ± 0.03) GeV 

C±{0) = 0.266 ±0.032 ^ 
fB = (192.8 ± 9.9) MeV 
fl = (220 ± 5) MeV 

«2 



0.08 ± 0.06 24 



24 



21 



17 



23 



14 



24 



m,{m,) = (1.27l°:°I) GeV 

C||(0) = 0.118 ±0.008 
fit. = (163 ± 8) MeV 
a^'" = 0.03 ± 0.03 1^ 
\b = (0.51 ± 0.12) GeV 



22 



17 



17 



Table 1 : Inputs used in the numerical analysis of the forward-backward asymmetry. All other 
inputs are taken from the PDG (22]. 



where Og is the angle between the £~^ and the B in the dilepton centre of mass system. Complete 
expressions for the differential decay width and the twelve /f functions can be found in Ref . y^ . 
The forward-backward asymmetry is obtained integrating Eq. (J9]) over cos 6f 
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(13) 



Using Eqs. ( [4||5|l2|13[ ) and the leading order expressions for the form factors, one finds that the 
numerator of the forward-backward asymmetry is proportional to 



mbniB 



eiRe C*, Cfiq') + 2'-^^^Cr + Oia,, A/m,) . 



q^ 



(14) 



At order 0{a^) this quantity does not depend on the form factor (^y. For this reason hadronic 
uncertainties affect the location of zero of the ^fb spectrum 
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(15) 



only at the NLO level. Note that this statement does not extend to integrals of the forward- 



backward asymmetry over an arbitrary q^ interval. In fact the denominator of Eq. (13), i.e. the 



differential rate, depends at leading order on both soft form factors through the trans versity 
amplitude A^l^r. at this order the form factors V ^ A^^ Ti and T2 are proportional to ^± and 
the form factor A2 is proportional to ^± — ^y. The uncertainty on our knowledge of these 
two form factors is the dominant source of uncertainty on the forward-backward asymmetry 
spectrum and various integrals. In Table [T] we summarize the inputs that we need for a complete 
numerical analysis of the forward-backward asymmetry; note that the C±,||(0) are taken from 
Ref. |14|p^ from where the q^ dependence of the form factors is also taken. 

In Fig.jljwe plot the differential forward-backward asymmetry in the 1 GeV^ < q^ < 6 GeV^ 
range. The blue shaded area in the first panel indicates the total uncertainty that we obtain 
by varying all the input parameters given in table [T] and the factorization scale /i^. The latter 
dependence is an artifact and can be reduced by including higher order terms of the perturbative 
as expansion. In this paper we allow /X5 to vary in the usual [m^^^/2, 2m^^^] range (but also 
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Figure 1: Error analysis of the forward-backward asymmetry spectrum. In the first panel we 
present the total error on the FBA spectrum obtained by adding all sources of uncertainty in 
quadrature. In the remaining panels we show the errors induced separately by the two soft 
form factors and by the variation of the low-scale /x^. 



present the results we obtain for the more conservative [4, 5.6] GeV range adopted in Ref. [14]). 
The location of the zero of the forward-backward asymmetry spectrum is c^ = (4.0±0.13) GeV^. 
In the other three panels we show separately the impact of the three most important sources 
of uncertainties. Note that the variation of the scale /x^ at which the Wilson coefficients are 
evaluated induces a constant shift on the whole spectrum and that the uncertainties on the 
g^ = values of the soft form factors switch sign above and below c^. 

The latter observation leads us to propose a new integrated observable for which we ex- 
pect a very tiny form factor uncertainty, namely the ratio between the the forward-backward 
asymmetry integrated in the [4,6] GeV^ and [1,4] GeV^ bins: 
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The SM predictions for the integrated forwards-backward asymmetry in the two low-g^ bins 
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Figure 2: Error analysis of the forward-backward asymmetry spectrum for C7 > 0. 

and for the ratio 7^[(?o] are: 

^fb(<?') dg2 = -0.089 (1 ± 0.085^, ± 0.17^^(0) ± 0.10^||(o) ± 0.031,est) 
= -0.089 ± 0.020 [22%] , 

^FB(<?')dg2 = 0.062 (1± 0.089^, ± 0.18^^(0) ±0.105„(o)±0.067,est) 

= 0.062 ± 0.015 [24%] , 

-0.701 (1 ± 0.19^, ± 0.012^^(0) ± 0.0057^1,(0) ± 0.083rest) 
-0.701 ±0.15 [21%] . 



^[4] 



(17) 

(18) 

(19) 

(20) 
(21) 
(22) 



where the suffix "rest" stands for variations of m^, rric^ /b, /x* ^ ^12 ^^d A^. Note, in par- 
ticular, that form factors uncertainties on TZ[q^] are vanishingly smaU while the impact of all 
other sources of uncertainty is enhanced. The basic reason behind this behavior is that the 
effect of any parameter with the exception of the two soft form factors is a parallel shift of 
the whole forward-backward asymmetry spectrum, thus enhancing their impact on the ratio 
Tl[ql]. Adopting the narrower variation /X5 G [4,5.6] GeV suggested in Ref. [14], the scale 
uncertainties in Eqs. (flTl [191 [2TI) shift from (8.5%, 8.9%, 19%) to (2%, 2%, 4%); consequently 



the total uncertainties in Eqs. ([18| [20| [22|) shift from (0.020,0.015,0.15) to (0.018,0.014,0.07). 
The impact of this more restricted scale variation is fairly small on the integrated FBA in the 
two bins because the total uncertainty is dominated by the form factors. On the other hand, 
57?.[4] is controlled by the scale variation and, as we explained above, the scale uncertainties in 



Eqs. (17) and (19) are anticorrelated. There are three main arguments in favor of adopting the 



observable that we propose. 

• Tl[q^] is sensitive to the same combinations of Wilson coefficients that determine the position 
of the zero of the spectrum but it is much simpler to determine experimentally. Moreover, 
it can be measured even in presence of new physics that removes entirely the zero from the 
spectrum. 
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The dominant source of uncertainty on TZ[q^] is the scale at which the Wilson coefficients 
are evaluated. Within the QCD factorization approach this dependence can be reduced by 
NNLO perturbative calculations and can, therefore, be brought under control. 

The branching ratio and the location of zero of the FBA spectrum are completely insensitive 
to the phase of the Wilson coefficient Ciq. In fact, the FBA is proportional to (7io and 
the location of the zero depends only on |Cio|. The integrated forward-backward asymmetry 
dependence on this phase is, on the other end, non trivial. This can be understood by looking 
once more at the LO expression for the numerator of the forward-backward asymmetry for 
arbitrary complex Wilson coefficients: 

^FB(g') oc Re 1^1*0 ('c9 + y(g^) + 2^^CTH (23) 

= Re jlCiol e-^^- (\C,\ e^^^ + Y{q') + 2^^|C7| e^^^) \ (24) 

= |Cio| cos((/)io - 09) IC'qI + COS 010 Re Y{q'^) + cos((/)io - ^7) 2 ^l^rl 

L q^ 

+ sin(/)ioImy(g2)l, (25) 



where we remind the reader that in the SM we have ^lo = (f)g = and 07 = tt. Eq. (25) 
displays explicitly the dependence of the integrated FBA on 0io. Unfortunately an effect 
here can always be reabsorbed into a change of the absolute value and sign of this coefficient 
(note that the FBA is actually proportional to |Cio|). We note also that the FBA is a 
parity odd observable that can be decomposed into CP even and odd components, the latter 



being the term proportional to sin0io. From the inspection of Eq. (25) we immediately see 
that the ratio TZlq^] is completely insensitive to \Cio\ (there is a mild residual dependence of 
TZlql] on |Cio| through the denominator of the FB asymmetry); in contrast, the dependence 
on the phase of Cio is quite strong and it arises through the interference with the strong 
phase in the matrix elements of the QCD penguin operators (the imaginary part of Y{q'^) 
for q^ < 4:ml originates from the QCD penguin matrix elements involving a light quark loop 
and is independent of q^ in the limit of vanishing light quark masses |25|). 

Let us consider the simpler scenario of vanishing new physics weak phases in C7 and Cg. In 
this case we have: 

AMq^) oc |Cio|cos0io[|C9|+Rey(g2) + 2^^^|C7|+tan0ioImy(g2)l . (26) 

L q^ J 

Integrating over the [1,4] GeV^ and [4,6] GeV^ ranges and assuming the SM values for the 
coefficients C^^q we obtain (formulas valid at NLO will be presented in the next section): 

^r,i Xf^FB(9')dg^ 1 + 0.13 tan 010 

7v 4 oc -^. oc — -— ■ — . (27) 

/'^FB(<?2)dg2 1-0.03 tan 010 ^ ^ 
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qr 


23.98 


-21.65+8.820 i 


35.42 


23.23+3.561 i 


qs 


1.159+1.260 i 


-1.120-0.1590 i 


1.798+1.859 i 


1.142+1.521 i 
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2.694 


9.969+0.3029 i 


1.993 


9.110+0.3954 i 
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-0.6372-0.2986 i 


-16.38 


1.294-0.3303 i 


-12.92 
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3.878-4.161 i 
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17.90 





12.03 


9710 


-5.718 





-8.446 





<?89 





1.042+1.088 i 





0.6379+0.6671 i 


<?810 


-0.2763-0.3005 i 





-0.4288-0.4432 i 





9910 


-0.6424 





-0.4754 






Table 2: Coefficients qi and qij. 



A very interesting feature of Eq. (27) is that, in this scenario, the difference between 7^[4] and 
its SM expectation is proportional to tan^ig. This behavior is usually displayed only by CP 
violating quantities; for instance, the CP asymmetry in the forward^backward asymmetry 
studied in Ref. |26| is proportional to Im Cio = |(7io| sin^io (see Eq. (6.10) in Ref. [26j). 
The ratio 7^ [4] retains this dependence on ^lo but is independent of |Cio|. Therefore we 
conclude that 7^ [4] not only complements the branching ratio and the location of the zero of 
the spectrum but also provide a unique and direct access to the phase of Ciq. 

Let us now discuss scenarios in which the spectrum of the forward-backward asymmetry does 
not contain a zero. In this case all sources of uncertainties simply shift the spectrum (see Fig. [2]) 
and the ratio 7^[4] turns out to be extremely clean. For illustration, we choose SC^dnt) c^ 1.15 
(this value reproduces the central value of the experimental determination oi B ^ -^sT) and 
obtain: 



(^[4]) 



Cr>0 



1.42 + 0.04 



(28) 



showing significant reduction in the error. 

Finally we summarize the present experimental determination of the B -^ K*i^i~ forward- 
backward asymmetry and attempt a first comparison with the SM prediction. The main limita- 
tion we encounter is that the Belle [12 j and CDF [l3j collaboration^ present results in the whole 
low dilepton invariant mass region (0 < mu < 8.68 GeV^) including regions (i.e. mu < 1 GeV^ 



''"The BaBar collaboration published a study of the integrated FBA in the whole low q^ region |11 , without 
giving separate results for the q^ bins we are interested in in this work. 



and mil > 6 GeV^) in which the theory approach we are adopting is subject to larger uncertain- 
ties. The Belle and CDF results, their weighted average and the corresponding SM predictions, 
with the caveat we just mentioned, read (the suffixes refer to the q^ bin expressed in GeV^): 





.[0,2] .[2,4.3] .[4.3,8.68] 
■^FB "^FB -^FB 


Belle 
CDF 

Average 


0.471^1 ±0.03 0.1llgi^±0.07 0.45lgi^ ± 0.15 

0.13lJ;^^±0.25 0.191°:^° ±0.14 -0.061°;^° ± 0.05 

0.45 ±0.28 0.12 ±0.27 0.24 ± 0.23 


SM 


-0.027 ±0.004 -0.050 ±0.012 0.133 ±0.026 



Given our ignorance of the correlations between the experimental uncertainties in the various 
bins, we are unable to extract meaningful constraints on the ratios we are interested in. 

3 Model independent analysis 

Let us start by presenting a numerical formula that allows to calculate the integrated observable 
i?[4 GeV^] for arbitrary values of the complex Wilson coefficients 07^8,9,10(^5) = Cff^i^^iJii)) + 
5C'7,8,9,io. We write: 



i?[4] 



/hn / / /ln 



10 



/hd V /ld . 

10 

A = 90 + 5Z Re[g, 5C,] + Y, Re[q,j 5C, SC] 

i<j=7 



(29) 
(30) 



i=7 



where in the suffixes AB, A = L,H stands for the low ([l,4]GeV^) or high ([4,6]GeV^) bin 
and B = N, D stands for the numerator or denominator of the forward-backward asymmetry; 
therefore /hn (/ln) and /hd (/ld) correspond to the integral of the numerator and denomi- 
nator of the forward-backward asymmetry in the [4, 6]GeV^ ([l,4]GeV^) range. The complex 
coefficients Qi and qij are given in Table [2[ 

In the following we will adopt the numerical formulas for BR{B -^ Xg^) and BR{B -^ 
Xs£^£~) presented in Refs. fSpO); the experimental determinations of these two branching ratios 
can be found in Refs. 27-34] and their world averages are given, for instance, in Refs. (5]- 



As we discussed in the previous section, although TZIqq] is an intrinsically CP conserving 
observable, the interference between weak and strong phases in the numerator of the forward- 
backward asymmetry introduces a strong sensitivity to new physics CP violating phaseq^ Espe- 
cially interesting is the sensitivity to the elusive phase of Ciq. We present the result of a model 
independent study of new physics contributions to Cr, Cq and Cio in Figs. |3}(4} When allowing 
for complex contributions to 5Ci^ we force the phases to lie in the [— 7r/2, 7r/2] range and write 
SCi = ±\5Ci\e^^' where (pi = arg^C^ or {diig SCi — tt) mod27r- In the upper panels of Fig. p^ we 



■'-The impact of CP violating phases in the denominator of the FBA, i.e. the branching ratio, is much smaller 
and essentially confined to the interference between C7 and Cg. 
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Figure 3: Upper panels. 95% C.L. regions allowed hj B ^ X57 in the [5C^^ 5C^] (for real 5Ci) 
and complex 5C^ planes. Lower panels. Dependence of 7^[4] on the Wilson coefficients C^. 



show bounds implied by 5 ^ X57 on the NP contributions SC^dut) and SCsif^b)- In the left 
panel we consider real coefficients and in the right one we entertain the complex SC7 scenario. 
The two distinct regions that are allowed in the left panel correspond to the two possible signs 
for the 6 ^ 57 amplitude that is approximatively proportional to C^dub) = C^^dub) -\-5C^{ii}))^ 
with C^^/Lib) < 0. The precise experimental determination of BR{B -^ Xg^) implies that 6C7 
can lie either in a small region around (same sign solution, C7 ^ C^'^ < 0) or in a narrow 
band around SC^ ^ 1 (opposite sign solution, C7 ^ —C^^ > 0). In the two lower panels of 
Fig. [3] we show the dependence of 7^[4] on the complex coefficient 5C7] the parts of the curves 
that are actually allowed by the B -^ X57 constraint are indicated by thicker lines. From the 
inspection of these plots we see that these effects do not depend much on possible NP phases 
in SC7 and that for the same sign solution C7 < 0, moderately large effects are possible. The 
opposite sign solution C7 > yields a change in the sign of 7^ [4] but is currently disfavored by 
B -^ Xsi^i- data. 

In Fig. |4] we present a similar analysis for NP effects on Cq and Ciq. Note that, in the 
absence of extra phases, 7?. [4] is almost completely insensitive to |5Cio|. Note the contrast 
between complex effects in Cq and Ciq. In the former case the impact of a non- vanishing ^g 
can be sizable but not essential in order to obtain large effects on 7^[4]; in the latter case, the 
condition ^lo 7^ is necessary in order to have an effect at all. 

An interesting feature of these plots is the appearance of very large effects on 7^[4] for certain 
values of the coefficients C^^q. Positive new physics contributions to C^^q tend to shift the zero 
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Rm 



15Ciol=4 



|d^Ciol=8/ 




^C9(//o) 



±I^Ciol 



Figure 4: Upper panels. Dependence of 7^ [4] on the Wilson coefficients Cg and Ciq. Lower 
panels. Bounds on the Wilson coefficients C9^io(/X5) from B -^ Xgii (dark green regions) and 
from a possible measurement of the ratio i?[4] (light green regions). We fix C7 to its SM value 
and assume i?[4]exp = ^[4]sm ± 50%. 
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of the FB asymmetry towards smaller q^ values. The integral of ^FB(g^) in the [1,4] GeV^ 
range receives, therefore, positive and negative contributions and, for some value of the Wilson 
coefficients, it can vanish. At this point the ratio 1Z[^ can become very large and one may need 
to adopt a different strategy (e.g. consider lZ[q^] for some other q^ value). 

It is particularly interesting to extract the bounds on 5Cq and SCi^) that follow from a 
hypothetical measurement of 7^[4]. In the lower panels of Fig. [4] we assume an experimental 
determination of the ratio 7?. [4] with an uncertainty of 50% and central values agreeing with 
the SM and study its impact on the real [5C9,5Cio] and in the complex 5(7io planes. In these 
plots the shaded disks represent the area allowed hy B{B ^ Xgii) and the light shaded areas 
are the region of the [Cg, Cio] plane that is selected by the measurement of 7^ [4]. 

From the inspection of the lower-left panel of Fig. |4] we conclude that even a relatively poor 
determination of 7^[4] can help hugely in determining the value of Cq. The lower-right panel of 
Fig. [4] illustrates the impact of a complex phase in SCiq: as long as ^lo 7^ a determination of 
7^ [4] strongly constraints |5Cio|. 

4 Few concrete examples 

In this section we briefly entertain three distinct new physics scenarios and investigate the role 



of a measurement of 7^[4] on their phenomenology. In Sec. 4.1 we consider the addition of a 
fourth generation to the SM. In Sec. |4.2| we discuss one supersymmetric extension of the SM: 
an i?-parity conserving MSSM with extra sources of flavor changing interactions in the down 



squark sector. In Sec. 4.3 we consider a family-dependent U{iy model. Note that in the 4th 
generation and Z^ models the effects are driven by new CP violating phases in Cio] in the SUSY 
scenario, the effect is controlled by more traditional contributions to the (chromo-) magnetic 
moment coefficients C^^s- 

4.1 4th generation 

The inclusion of a fourth generation is perhaps the simplest extension of the Standard Model. 
The phenomenology of this model, usually referred as SM4, has recently been the subject of 
intensive scrutiny (see, for instance, Refs. |35j - |49| ). In this work we are interested in possible 4th 
generation effects on B ^ K*£i^ speciflcally on the ratio 7^[4]. The two parameters that control 
SM4 effects on the magnetic and semileptonic coefficients are the mass of the 4th generation 
top quark, m^/, and the ratio of CKM4 elements A^^, = V^/^V^/^/V^^V^s- From the analyses in 



Refs. 44,45] we extract ranges for these two parameters that are phenomenologically viable; 
in the numerical analysis we take m^/ = 500 GeV and |A^^/| < 0.1. 

The SM4 expressions for the f contributions to the Wilson coefficients C7-10 can be ob- 
tained trivially from the corresponding SM expressions. Since we are interested in the forward- 
backward asymmetry in B ^ K^'M^ we do not need separately the combinations V^/^V^/^ and 
V^^V^*, but only their ratio. 

In Fig. [5] we show the dependence of TZ[^] on the magnitude and phase of A^^, for m^/ = 



500 GeV. It is clear that large deviations are possible though only in presence of large phases 
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arg At't.= 180 
arg AJ,= 
arg 4'- -120 



R[4] -oi 
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-— |A,^,h0.01 



Figure 5: Expectations for 7^[4] in SM4. We fix the mass of the t' to 500 GeV and aUow the 



ratio of CKM4 elements A^^, to vary according to the bounds obtained in Refs. 144,45 



in the new CKM4 elements. 

This strong dependence on the CP violating phase of A^^, stems entirely from the 4th gen- 
eration contributions to Cio- In fact, the functional dependence of SCf'^ and SCf^ on rrit' is 
such that for large t' masses (in the range [200, 1000]GeV) one has 2 < \5Cf^/5Cf'^\ < 10. 



4.2 MSSM 

Let us consider an i?^parity conserving MSSM with non- vanishing sources of fiavor changing 
neutral currents in the down squark sector (i.e. non Minimal Flavor Violating). As an example, 
we focus on contributions induced by a non-vanishing value of the LR mixing between the 
second and third generation. Adopting the formalism of Ref. 50 we define this mixing in the 



super-CKM basis (a basis in which the squark mass matrices are subject to the same rotations 
that diagonalize the regular quark sector). 

In this model, gluino vertices are fiavor changing and the dominant contributions to the 
Wilson coefiicients come from 1-loop involving gluino and down squarks. The parameter that 
we need are the gluino and down squark masses and the mass insertion (62^) lr (see Ref. 50 
for more details). 

The dominant constraints on this insertion come from B -^ -^s7, B -^ Xsi£^ Bs mixing, 
vacuum stability and absence of color breaking minima in the potential. The interplay of these 
constraints and the allowed ranges for {S2^)lr have been studied at length in the literature. 
In this work we use the results of Ref. [5l]. Note that in Ref. 51 the bounds are given for 
rriq c^ rrig c^ 350 GeV and can be easily rescaled to cover the rriq c^ rrig c^ 500 GeV that we 
consider here. 

In Fig. [g] we show the size of supersymmetric corrections to 7^[4]. The opposite sign solution 
for the B -^ Xg^ amplitude (i.e. the C7 > scenario) requires a very large mass insertion, 
(^23)^-^^ ^ —0.08 X {rriq/ 500 GeV), that is excluded by the charge and color breaking bound 
that reads \{S2^)lr\ < 0.017 x {rriq/ 500 GeV). On the other hand, the dependence of 7^[4] on 
the phase of the mass insertion is not as strong as in the 4th generation case. The reason for 
this is that in this supersymmetric model we do not have large contributions to Cio] therefore. 
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Figure 6: Dependence of 7^ [4] on |(^|^^| for rrig c^ rriq c^ 500 GeV. We indicate with a thicker 
hne the region aUowed hy B ^ Xgj. Sohd, dashed and dot-dashed hnes correposnd to various 
values of ^^^ and 



the mechanism (see Eq. (27)) that induces a large sensitivity to CP violating phases in 7?.[4] is 
not active in this case. 



4.3 Family-dependent U{iy 

As a final example of new physics scenarios in which the observable 7^ [4] has a strong dis- 
criminating power, we consider a Z^ model with tree-level fiavor changing couplings proposed 
and analyzed in Ref. 52 (see also Ref. 53 for an analysis of an eflFective fiavor changing Zbs 
coupling). In this analysis we assume that the U{iy coupling constant {g2 in the notation of 
Ref. 1 52 1) is identical to the SM U(l) coupling gi and that the Z and Z^ have identical couplings 
to leptons. Therefore, the only parameters that remain (for what concerns b ^ s transitions) 
are the Z' mass and the fiavor changing complex couplings B^^ , defined as 



Cfc = -g2Z'^h'[BtPL + B^^Pn]s 



(31) 



The efiFect of 5^ is to generate contributions to the operators Oq ^q obtained from Og^io by 
replacing the left-handed quark current with a right-handed one. In our study we will not con- 
sider this possibility and assume that Z^ couplings follow the weak interactions chiral structure 
(i.e. we take 5^ = 0). Utilizing Eq. (31) we obtain: 



SCio 



4^2 



1 



%,%: 



.Ctp 



4:71 



92mz f.L 
giTUz' 



1 
2 



VtbVu 



.OLf. 



47r 



g2'mz 
gimz' 



B 



bs ") 



(32) 
(33) 



where sw is the weak mixing angle and, following Ref. |^2|, we defineqjS^ = B^^{g2mz / QiTu'^) . 
The reason for replacing S£ with the rriz' dependent coupling 5^ is that the latter is directly 



^In Ref. 52 the rescaling of the B couplings is performed but the hat symbol is not explicitly introduced. 
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Figure 7: Expectations for 7^[4] in a family-dependent U{iy model. We fix the mass of the Z' 
to 1 TeV, allow the complex parameter |5£| < 2 x 10~^ [52] and assume that the Z' coupling 
to leptons is identical to the SM. 



constrained by data on B^ mixing and B -^ {n^ip^cp^rj^p^uj^ f^)Ks decays: the authors of 
Ref. 1 52 1 find that these constraints result into the upper bound \B^J ^ 2 x 10~^. From the 



inspection of Eqs. (32) and (33) it is clear that the exchange of a Z' generates only sizable 



contributions to Ciq. In fact, as a direct consequence of our assumptions of identical Z and Z' 
couplings to leptons we have SCq/SCiq = 4^^ — 1 c^ —0.075. For this reason we expect large 
contributions to 7^ [4] only for (large) non-zero values of the 5^ phase. 

We present the results of our numerical analysis in Fig. [Tj where we show the dependence of 
TZ[4] on the modulus and phase of S^ for m^/ = 1 TeV. Clearly very large effects are possible 
and a measurement of 7^[4], together with CP violation in Bg mixing and b -^ sqq modes can 
provide a smoking gun for these models. 



5 Conclusions 

In this paper we introduce a new observable in exclusive B -^ K''£^£~ and Bg -^ (j)£^£~ decays: 
the ratio 7^[4] of the integrated forward^backward asymmetry in the [4, 6]GeV^ and [1, 4]GeV^ 
bins. The separation between the two bins roughly coincides with the position of the zero of 
the FB asymmetry spectrum (calculated in the SM at NLO). 

We have shown that the bulk of the theoretical uncertainty on 7^ [4] is due to perturbative 
scale uncertainties and that the form factor dependence drops out almost completely. This ratio 



is therefore extremely interesting because its SM prediction Eq. (22) can be systematically im 



proved without requiring any non-perturbative input. The position of the zero of the spectrum 
shares these positive features and can be determined with even better theoretical accuracy, but 
is much harder to measure. 

A very interesting feature of 7?. [4] is the strong dependence on the CP violating phases 
(especially those appearing in the Wilson coefficient Cio). This dependence arises because 
of the presence of strong phases in the b -^ s£^£~ matrix elements of current-current and 
QCD penguin operators (i.e. Oi-q). The latter can be determined perturbatively using QCD 
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factorization (or SCET). 

We performed a model independent study of new physics contributions to the coefficients 
Cr^g^io and speciahzed the analysis to several concrete new physics models. We find that, after 
imposing the constraints from experimental data on radiative and semileptonic rare decays, 
there are two scenarios in which large effects are possible. 

• Models with modest (large) contributions to the coefficients C^ (Cg). As a representative of 
this class we considered an MSSM with non-vanishing ((^|^)Li? mass insertion. In this case, 
large effects on 7^ [4] and on the location of the zero are correlated. Moreover, we do not find 
a very strong dependence on the new physics phases of ^Cr^g. 

• Models with large and complex contributions to Ciq. We considered the inclusion of a se- 
quential 4th generation to the SM (SM4) and a U{iy model with family-dependent couplings 
(i.e. tree-level fiavor changing Z' interactions). In SM4 we find large contributions to Cio 
because of the loop-function dependence on the t' mass; in the family dependent U{iy model 
this behavior is a result of our assumption of universal Z and Z' couplings to leptons. In 
these models we find very large effects on 7^ [4] driven by the phase of 5(7io. At the same time 
there is almost no change in the position of the zero of the spectrum because new physics 
contributions to Cio do not impact this observable at all. 

In conclusion, we believe that this new observable should be promptly included in future exper- 
imental studies oi B ^ K'^ii decays because it allows for an easier access (i.e. less luminosity 
is required) to the new physics probed by the zero of the spectrum and, at the same time, is 
sensitive to the quite elusive new physics phase in Ciq. 
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